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Square unit in Husimi lattice (the Bethe square lattice) is generalized to be rhombus with ran¬ 
domly variable angles. The independence feature of unit cells in recursive lattice makes the random 
angle conformation possible in the model construction, which is unfeasible in the conventional lat¬ 
tice. Since the randomness of the conformations in real system is naturally introduced into the 
model, this new lattice methodology can describe the off-crystal metastable states without artificial 
randomness. With reasonable simplification, a coefficient A(0) is formulated to present the effect 
of angle in the rhombus unit. A “visit and count” recursive technique is developed to numerically 
calculate the thermodynamics. While the computation randomizes a quenched configuration in each 
iteration, the calculation counts and averages a large number of random units to deal with a system 
in equilibrium with annealed randomness at particular temperature. The critical temperature T c 
of spontaneous magnetization transition is lowered with the presence of angle randomness, which 
implies a less stable system. Besides consistent results to the regular lattice, the random-angled 
lattice features a distribution of solutions and the thermal fluctuation with exact calculation. The 
effects of the variation of energy and ground state parameters on thermodynamics are investigated. 


I. INTRODUCTION 

Except some well-arranged periodical conformations of particle in nature, e.g. the crystal structure, lattice models 
are usually employed in physics to quantize the continuous space for convenient and simplified modeling [lj]. With 
appropriate approximation, the lattice methodology works well to describe real systems or to investigate the physics 
insight, classical examples include the Flory-Huggins model to describe the thermodynamics of polymer solution^, 3|, 
Ising model [-Jill, Potts model Q and so on. Nevertheless the balance between the lattice approximation and the space 
continuity of real system is always an important concern in physics modeling. Among a number of sacrifices of lattice 
modeling, one notable compromise is the variable angle between interacting bonds. 

Imagine a free-arranged 2D interacting monatomic system, in which a particle interacts with four other neighbors, 
then a quadrilateral obtained by sampling any four neighbor particles has a random-angled conformation with variable 
energy, while in a regular square lattice any four neighbor particles must form a square with the angle 90°. Although 
the lattice methodology can count in the angle effect under some circumstances, for example to adjust the related 
bond interactions to count the effect of off-90° angle, or non-square quadrilateral with particular angles can serve as 
the basic unit to construct a lattice in a periodical fashion, the problem is the uniformity of lattice cells: however 
has the lattice been setup to fit the angle effects, the fitting must be a uniform property of all the unit cells, while 
ironically the angles actually exhibit the feature of randomness in real system. 

According to the above concern, we developed a lattice model to handle the random angles via the unique properties 
of recursive lattice. Recursive lattices such as the classical Bethe HI and Husimi [§, 101 lattice have been studied 
for decades and proven to be a reliable methodology to describe the real system or to approximate the regular 
lattice [III E3 ■ In a recursive lattice, the unit cells only interact to each other on the joint site and the conformation 
of one cell is independent, which makes individual conformation with variable angles possible. The ensemble of 
individually-shaped units raises a non-uniform lattice to better demonstrate the space randomness in real system. 
In this random lattice, the conformation of each unit is a rhombus randomized and fixed by the program, therefore 
locally it has a quenched randomness, nevertheless the calculation counts and averages the statistical variables over 
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a large number of random units to deal with the system in equilibrium at particular temperature (details will be 
discussed below), subsequently the system can be treated as with annealed randomness in general. 

Ferromagnetic Ising model on this new lattice is chosen to investigate the thermodynamics. The first reason is that, 
it is a classical system that has been throughly studied for many years fT?fl42lT . we can easily compare and validate the 
calculation to previous studies. Secondly, this is an exactly solvable thermodynamic model [la I2ai32] . and the thermal 
fluctuation indicated by the exact calculation imprints the point of stochasticity in an exactly-solvable deterministic 
system that we would like to address in this work, i.e. the thermodynamic calculation still follows a deterministic 
way however it gives out the randomness in the real system. And another significance of applying Ising spins on this 
random lattice, is that we can combine the idea of spin system and structural glass model to investigate the off-crystal 
metastable state. Since the amorphousness of glass can be presented by the random conformation of lattice, the spins 
can form a glassy state with their well-defined orientation, that is, a glass of spins without the concept o f “sp in glass”. 
This model is then able to describe some particular systems, for example the metallic glass ferromagnet |33|, in which 
the materials is in glassy state but the well-defined orientation of atoms gives superparamagnetism. 

Other than the application in metallic glass ferromagnetism, we believe this methodology could be a decent contri¬ 
bution to this field, and has the potential to be applied in various physics problems, for example the glass thermody¬ 
namics, vibration of crystal lattice, and quasicrystal. In this paper we will explore the random-angled conformation 
in Husimi square lattice. More specifically, the unit cell is generalized to be rhombus with only one independent angle 
variable for convenience. It should be easy to extend our work to be a more general quadrilateral lattice, or other 
polygon and 3D cases. The Hamiltonian formulation and thermodynamic calculations on this lattice will be discussed 
in details, and we will investigate the thermal behavior, phase transitions and the effects of energy and ground state 
parameters setup. 


II. LATTICE CONSTRUCTION 

The conformation of a unit cell is independent in a recursive lattice, therefore a Husimi square lattice can easily 
transform to be an recursive ensemble of random-angled rhombus units. FigJTJi shows the idea of the recursive rhombus 
lattice, the structure is still in a Husimi recursive fashion but with the basic unit cells as random rhombuses. In a 
more general view, the basic unit cell of this recursive lattice can be called as a “conformation distribution cloud” 
instead of a fixed geometry unit (FigO})). The “cloud” collapses to be a particular rhombus conformation when it is 
visited during the calculation. 

Note that in the lattice all the edges are in the same length to represent the neighbor interactions, the different 
size in FigQJi is for possible drawing on the plane page. With the fixed edge length only one independent angle 9 is 
necessary to determine the structure of a rhombus. By defining 9l as the lower limit, we have 9 £ [9l : 180° — 9^\. In 
this paper, the angles are randomized with even probability distribution in the region [ 9 l, 180° — 9l\, and 9l is set to 
be 20°. 

For the angle 9 between two bonds in the rhombus, we must have one or more the lowest energy angles denoted as 
9 0 at which the E(0 O ) = E m i n , i.e. the most stable state. For example with 9 a = 90° the square unit will be the most 
preferred configuration. In this work we will only discuss the case that 9 0 has a single value in the range (0, 180°). 
However a multi-stable states system implied by many values of 9 a is possible. 


A. Modeling Angle Energy 

With the ground state angle 9 0l a higher energy E(9) of the angle 9 £ [9l,tt — 9l\ rises with the deviation from 9 a , 
and then we should have E(9) = f(9, 9 0 ). 

Previous studies on Husimi lattice usually consider only the interaction J between nearest-neighbor sites, while 
in several works the model also counted in other interactions, e.g. the diagonal interaction j3j|. Following this 
consideration, by safely assuming that the angle does not affect the nearest-neighbor interactions, then if the angles 
would do anything on the system energy it should be either the angle bending energy itself, or effects on the diagonal 
interactions in the rhombus, which relates to the opposite angles with the law of cosines. Since an artificial bending 
energy on the angle is similar to an external field applying on the lattice sites and not quite interesting, we will take 
the second choice to formulate the angle energy. 

In homogeneous lattice the interaction energy of a pair of diagonal spins ( Sj , Sj>) is E p = — J p ■ SjSj’, where J p 
is the diagonal energy parameter. Assume the diagonal interaction is linearly proportional to the diagonal length, 
according to the law of cosines c 2 = a 2 + 6 2 — 2 abcos9 1 with c as the diagonal opposite to angle 9 and the side length 
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FIG. 1: (a) A sample demonstration of Husimi rhombus lattice with random-angled units, (b) The basic rhombus 
unit with random angle in the range of [Ol, 180° — 6l\ with the lowest energy angle 9 0 , in this work the 6l is set to 

be 20° unless specially notified. 


of rhombus a and b to be unit 1, we have 


E p (0) = -J p 


SjSj' 


^/2(1 — cos0) 
\/2(l - cos 6 0 ) 


Define A(6) as the angle energy coefficient: 


( 1 ) 


m = 


1 — cos 9 

1 — COS do ’ 


( 2 ) 


as a natural choice for the quadrilateral, in this work we take 9 a = tt/2 for the general discussion. In this way, we 
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then have 


E p = -J p -A{0)-S j S r . (3) 

where A{9) = \/l — cos 9 for a random angle 9. The effect of various 9 0 will be discussed later. 

It should be addressed here that the above formulation of angle energy is counter-intuitive because the assumption 
“diagonal interaction linearly proportional to the length” implies that the farther the diagonal pairs particles are, 
the larger the interaction is, which is not a common case in nature. The reason of such formulation is that when we 
combine the energies of two diagonal pairs in one rhombus, a higher total energy (smaller weights) of any off-angle 
configuration can be guaranteed with a larger magnitude of energy change on the other term: J p ■ A(n — 9) . In another 
word, the formulation based on A oc yj [1 — cos 9 a )/(I — cos 9) will lead any off-0 o configuration corresponding to a 
larger energy (lower weights), and the structure with 9 = 0 O is the most stable configuration with the largest weights. 

We formulate the angle energy in this way merely for convenience to preliminarily explore the nature of random- 
angle lattice. It should be easy to extend the methodology with other angle energy formulations for particular 
circumstances. 


B. The Hamiltonian 


Adding the angle term introduced in Eq|3] the energy of a lattice unit containing four magnetic energies, four 
neighbor and two diagonal interactions is: 


e 


4 

- E J - s ^ 

<i,j> 




J p (A(9i) ■ SjSj. + A{0j) • SiSi+2), 


( 4 ) 


and its Boltzmann weights is, 


w ( 7 ) = e /3e , 


( 5 ) 


where 7 denotes the conformation of unit. 

In this paper the interaction parameter J will be set as +1 for the ferromagnetic case. Being consistent to J, J v 
is always set as positive. The effect of J p variation from 0.1 to 0.6 will be discussed and J p = 0.2 is taken to be the 
reference case. The magnetic field is not concerned with H = 0. In the third term Sj and Sji denotes to the spins 
neighbor to Si and they assemble a diagonal pair with opposition to the angle 9i on the site Si . The Hamiltonian of 
the entire lattice is the summation over all the units energies 

E = 'Y^e a , ( 6 ) 


where a denotes the index of unit cells. 


III. CALCULATION 

The calculation principally follows the same methodology of our previous works [29l [3lL l32i. 134] on Ising model on 
Husimi lattice. Imagine the entire lattice has an original site So on the 0th level, then the sites nearest to So are 
marked as Si (or S'i to distinguish the two sites in the same unit), the site diagonal to So is marked as S 2 , and so 
on. The situation of ?zth level is shown in the Fig(2] On one rhombus unit there are three sub-trees Tn+i and T~ n +2 
contributing to the sites S n , S' n +i and S n + 2 , and a larger sub-tree T n is synthesized by adding the weights of the 
local unit. With the particular spin state of S n = ±1, define the partial partition function (PPF) Z n (S n ) on level n 
to count the contribution of sub-tree to the site S„ (not included), then the nth level PPF is the function of PPFs on 
the previous level n + 1 , n + 2 and the weights of local unit u;( 7 ): 


8 

Z n (+) = y^ j Z n+ i(S n+ i)Z n+ i(S , n+1 )Z n+ 2(S n+ 2)w('y), 

7=1 


( 7 ) 
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FIG. 2: Four sites Si in a rhombus unit are labeled with levels. The 71 denotes the sub-tree contributions from level 
i to the local unit. Solutions Xi on corresponding site Si are calculated accordingly. 


16 

Zn(~) = J2 Z n+l{S n+1 )Z n+1 (S' n+1 )Z n+2 (S n+2 )w( 7 ). 

7=9 

And the partition function (PF) of the whole sub-tree on nth level is given by 

= Z n (+)e pHS “ + Z n (-)e~ pHS \ 

which summers all the weights at level n including the spin S n . 

By introducing the ratios 


fin 


z n (+) 

Z n (+) + Z n (—) ’ 


( 8 ) 

(9) 


( 10 ) 


and a compact note 


Vn = 


Zn{ ) 

Z n {+) + Z n {—) 1 


i{S n ) 


X n if S n = +1 
y n if S n = 1 


( 11 ) 


( 12 ) 


As the weights ratio of PPFs with the spin state of S n , x or y denotes the sub-tree contribution to the magnetization 
of the site S n , i.e. the cavity contribution. Therefore, we can use these ratios, which we can the “solution” of the 
model in the following, to indicate the magnetization of a site. In terms of 


B n — Z n (+) + Z n (—), 

we have Z n (+) = B n x n and Z n {—) = B n y n . With Eci 171 l8l and fl2l it gives 


B n ^ni^S) — ^ ^ B n +iZ n -\-i (<Sn-}-l)-Bn+l~re-|-l (S n+1 )B n+2 z n+2 (>S f n _(-2)'ic(y), 
*.(±) = E z n+ i(S n+ i)z n+ i(S' n+1 )z n+2 (S n+2 )w('y)/Q(x n+1 ,x n+2 ), 
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where the sum is over 7 = 1 , 2 ,3,..., 8 for S n = +1, and over 7 = 9,10,11,..., 16 for S n = —1, and where 

Q(^n+l 7 ^n+2) = [-^n/-^ n -|_i-^n+ 2 ] 5 


it relates to the polynomials 


8 

Q-\-{%n-\-l i ^n+2) = ^ ^ -£72+ 1 ( *$ti +1 ) • 2 - 71 +l(<Syx-|-i)-£71+2(^72+2) ^(7)7 
7=1 
16 

(^m+1 7 ^m+ 2 ) = ^ ^ n+l(^n+l)^-n+l(<S > 77,-|-i)^-n+2(*^ > n+2)'^(7)7 
7=9 

according to 


Q(*£n +1 7 *£71+2) — Q+ (^n +1 5 *£72+2) Q — (^/ti+I 5 2/71+2)- 

In terms of the above polynomials, we can have the recursive relation for the ratio x n in terms of x n +i and x n+ 2 : 


^ _ Q-\- (*^ 71+1 1 *£71+2) 

^(^ 71+1 5 ^n+2) 

Similarly for y n = 1 — x n : 

_ Q— (2/71+15 2/71+2) 

Vn Q{yn+l,yn+2) 

Therefore the solution x is a function of PPFs and subsequently a function of zs on the previous levels: 

*" = z n (+) n i + z n (-) = f{Xn+1 ' Xn+2) (15) 

then starting with two initial guesses of xs we can recursively calculate the x on lower level. 

In the homogeneous lattice the recursive calculation is done for many iterations until reaching a fixed cycling solution 
%k = Xk+i (the fixed point solution). However in the random lattice discussed here, the feature of configuration 
randomness and variable local weights uj( 7 ) in each iteration makes the exact fix-point solution unachievable. After 
a number of iterations the calculation will provide an oscillating solutions instead of the fixed x. Therefore the 
calculation is designed in a “visit and count” way: with a known solutions x n +2 and x n +i from previous levels, the 
calculation process “visits” the rhombus on the level n with base site S n , where the new rhombus is randomized out 
with particular angles when the “visit” occurs. (The index of levels is counted down because we assume the calculation 
is heading to the imaginary origin point, i.e. the level 0 of site So)- Assume the three contributions from sub-trees 
T n +i and T n +2 on previous level are identical j35|, the new solution x n on the site S n can be calculated by Eq[7] [ 8 ] 
and M with the random local weights + 7 ). After this “visit”, the solution obtained in this iteration is numerically 
averaged to the solutions on previous levels, and the x is the description of the entire lattice on level n. Subsequently 
we will have a sub-tree contribution T n for the next “visit”. Notice that although one site at level n is both on the 
levels visited and not yet visited, the visitation only determines the conformation of next random rhombus, the x n as 
a function of previous levels x n +i and x n +2 is not affected by the lower level, and indicates the probability occupied 
by a +1 spin on that site. 

According to our experience the calculation reaches the “stable” solutions oscillation x n £ ( x±cr) (a: the standard 
deviation) after no more than 2000 iterations. In the program we do the calculation 6000 times and average the last 
2000 ss to reach a reliable x. 

In this way although the exact calculation is executed in each iteration, instead of exact solution what we actually 
obtained is a solutions distribution. FigOJi shows a reference average solution with J p = 0.2 and 0 o = 90° in a wide 
temperature range, Fig [3jo shows a distribution of twenty sample solutions at T = 2. The temperature is normalized 
with Boltzmann constant to be 1. 

The “visit and count” method implies that the calculation process deals with a quenched randomness while a 


(13) 


(14) 
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FIG. 3: (a) The reference average solution x with J p = 0.2 and 9 0 = 90°. (b) The distribution of twenty sample 

solutions x at T = 2.0. 


local unit configuration is randomized out by program and fixed thereafter. Nevertheless the x averaging over a 
large number of lattice units describes an entire system in (or near) equilibrium, therefore the randomness of xs 
and subsequent thermodynamic properties can be treated as annealed. Otherwise it makes no sense to average the 
variables Xi on sites i for a finite pure quenched system, although the average of extensive properties such as free 
energy, which will be discussed below, is still feasible [38|. 


A. Calculation of Thermodynamics 

The Helmholtz free energy as the function of PF is given by F = —T\ogZ. Although the lattice is infinitely large 
and it has no sense to calculate Z or F, the free energy of the local unit on nth level can be achieved by the difference 
of free energies on successive levels Fiocai.n = F n (Z n ) - 2F n+1 (Z n+1 ) - F n+2 {Z n+2 ). 

Then we have 

Fiocai.n = -Tlog 2 Zn - 

^n+l^n+2 


(16) 
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FIG. 4: The thermodynamics of reference solution x with J p = 0.2, 9 0 = 90° 


With the derivation on PF, PPF and xs introduced above, the local free energy Fi oca i >ra can be simply calculated as a 
function of solutions x n , x n +i, x n +2 and temperature T. The same “visit and count” method is applied on the free 
energy. In an iteration once the F loca ] rl was obtained, it is averaged back with previous results: 


_ 1 F n+ 1 • (n - 1) + Fjocal.n 
2 n 


(17) 


The 1/2 is to divide the number of sites in one unit (four half-shared sites). Again in the calculation practice, we 
average the last 2000 results out of 6000 iterations in the program. 

The energy per site, i.e. energy density is calculated with the conformation probability. A rhombus unit with 4 
spins on each site has sixteen possible conformations 7 . The probability of one conformation is 

r>l \ \Zn+\(Sn+\)Zn+\{S' n +\)Z n +2(.S n J t -2)\w{pf) 

Fyy) = -~-; (,i») 

Then the energy density is defined as the summation over the product of the energy and the probability of a confor¬ 
mation state: 


E = ^E e (7)-P(7) 


where e is from the EqfU 

The average entropy per site is derived from F and E: 


S = 


E-F 

T 


(19) 


( 20 ) 


A reference thermodynamics of J v = 0.2 and 9 0 = 90° in a wide temperature range is shown in Figj4j The transition 
can be observed on the entropy or energy curve, on which a slope change at T = 3.096 indicates a phase transition 
from disordered to ordered state with decreasing T. When the T is close to zero, the entropy approaches to zero as 
expected and the energy or free energy at zero is given as ~ 2.180. This ground energy quantitatively measures the 
effect of angle/diagonal energies presence. Referring from our previous work[? ], the recursive lattice of coordination 
number q has the ground energy E{T = 0) = ^ J • q if there is no other interaction terms. Therefore the contribution 
of angle/diagonal energies to the ground energy is ~ 0.180. 
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FIG. 5: The comparison of the solutions around T c of stochastic and homogeneous lattice. The temperature 

resolution is AT = 0.0001. 


IV. RESULTS AND DISCUSSION 
A. Spontaneous Magnetization and Critique Temperature 

Fig[3ji shows clear spontaneous magnetization transition[U [37] with the presence of angle randomness. In the 
high temperature range the solution is 0.5 indicating that all sites have 50% probability to be occupied by either 
+1 or —1 spin (disordered state). With cooling process till the critique temperature T c , the solution decreases and 
approaches to zero, i.e. at low temperature the system prefers more -1 spins and in the near-zero range all the sites 
will be occupied by —1 spins (ordered state). It should be clarified that although in Fig[3Ji the transition is very 
sharp between T = 3.0 and T = 3.1, it is not a discontinuous transition and the detail is shown in Fig[5] Depending 
on the initial seeds of x the solution below T c may either symmetrically bend up and approaching to 1. That is the 
magnetization to the other direction and has the same thermodynamics therefore it is not shown in the figure. 

In the quantitative aspect, presence of angle randomness decreases the T c . FigOshows the solutions of deterministic 
case and one sample of stochastic case zoomed in around the transition range. The deterministic solution is calculated 
with the same formulation except setting 9 always to be 9 a in the program. The T c _^ etel is 3.131 and T c _ stoc h is 3.096. 
This is because 0 o corresponds to the lowest energy rhombus conformation, subsequently the homogeneous lattice 
with universal 9 = 0 o must be the most stable system. Meanwhile the presence of off-angle configurations, even the 
larger the off-angle degree the smaller weights they have, will generate higher energy to unstablize the system to be 
easier to experience the phase transition, i.e. the lower T c . 

Another interesting phenomenon is that the stochasticity does not randomize T c . Five samples of stochastic solutions 
in the around T c range are shown in Fig[6] While below T c all sets of solutions behaviors as a random distribution, 
the transition temperatures where solutions converge to the 0.5 are the same. This phenomenon can be understood in 
EqQ]and[5] the driving force of temperature applies on all the energy terms in Eq0]to obtain the weights, no matter 
what the randomized angles one rhombus unit have, the weights preferences on particular configurations under certain 
temperature are the same. In another word, regardless of the shape of the rhombus, when temperature decreases till 
the weights of the three and four —1 spins on the four sites are higher, the spontaneous magnetization occurs and the 
solution x heads to 0, or vice versa. 


B. Thermal Fluctuation 

From Figj5]we can more clearly observe that while the deterministic solution is a fix-valued line along the tempera¬ 
ture axis the stochastic case provides a distribution. The solutions fluctuates off the equilibrium state (the equilibrium 
state is an imaginary line assembled by the solutions x averaged over infinite xs , it can be expected to be a single 
line likewise the deterministic solution, however cannot be calculated and displayed in the figure). Furthermore the 
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FIG. 6: Five sample solutions in the around T c range. The temperature resolution is AT = 0.0001. 



FIG. 7: The entropy of J p = 0.2, 9 0 = 90° with resolution AT = 0.001, the entropy of homogeneous lattice is 

displayed for comparison. 


thermodynamics obtained from solutions also exhibit thermal fluctuation. More specifically it is possible to have a 
higher entropy at lower T, which is not possible in the exact calculation on homogeneous lattices. 

Unless displaying several samples together on the same graph, observing the fluctuation with large temperature 
increment is difficult, e.g. in the Figj3]and Figj4] In the following discussion we call the T increment as “resolution”, 
and usually we set AT < 0.01 to observe the fluctuation/distribution. The fluctuation degree can be characterized 
by a “middle-value-based” deviation which will be detailed in next section. 


1. Entropy Behavior 

Figm is the entropy with J p = 0.2, 9 a = 90° and resolution AT = 0.001, the entropy of homogeneous lattice is 
also displayed for comparison. Besides the similar behavior, entropies of two systems distinguish from each other as a 
distribution and fix-valued neat line respectively like the solutions in Fig[5] In the stochastic case we can locate many 
neighbor pairs of data points that the entropy value at lower T is larger than the other at higher T, this thermal 
fluctuation at equilibrium is not treated as disobeying the second law[36j and can be observed in both experiments 
and simulations. Nevertheless exact calculation usually cannot present the fluctuation and that raises an advantage 
of our model. 
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FIG. 8: Six samples of entropies exhibiting fluctuation and negative in the near zero range. 


2. Near-Zero Temperature Range 

The fluctuation exists rigorously in the range from a T slightly higher than zero to the T c . The solution is consistent 
0.5 above T c and 0 in the close-zero area, in these two ranges the driving force of temperature is too strong to allow large 
fluctuations. Within the fluctuation range, there is a more significant fluctuation found in the near-zero temperature 
range. Fig J8] displays six sample solutions with the resolution AT = 0.01, the calculation obtains negative entropy in 
a few higher T range. Even it is easy to explain this behavior by the randomness in the free energy and energy density 
calculation, this phenomenon is still interesting if we consider that the drive force of temperature rapidly pulls up the 
entropy away from zero above a bit higher T (T > 0.55), and firmly fix the entropy down to zero below a bit lower 
T (T < 0.20), but leave a range in which the randomness leads to a fluctuation with negative entropy. Several trials 
have been done and it seems that this negative S range is inevitable. The reason of this phenomenon is not clear yet. 


C. Effect of J p 


The effect of J p in homogeneous lattice has been discussed in reference [31]. Generally the J v of the same sign of 
J enhances the weights of same spins configuration and increase the stability of system. In the random lattice this 
principle also works. The T c variation with J p = 0.1 to 0.6 are summarized in table [i] 

TABLE I: The effects of J p on T c in homogeneous and random lattice. 


Jp 

^c deter 

-^c stoch 

0.1 

3.105 

2.949 

0.2 

3.132 

3.096 

0.3 

3.307 

3.261 

0.4 

3.475 

3.414 

0.5 

3.638 

3.564 

0.6 

3.799 

3.710 


The effect of J p on the fluctuation degree is a more intriguing finding. A larger J p will stress the weights of random 
terms in EqQ]thus we can expect a more rigorous fluctuation. For a single set solutions we developed a “middle-value- 
based deviation” S m to quantitatively characterize the fluctuation degree. For the solutions series {^(T))}, the middle 
value Xi_ m is calculated for each pair of two successive solutions x t and aq + i: Xi_ m = A+£i±±, then the deviation S m 
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FIG. 9: The solution distributions in around T c range of J p = 0.1 and J p = 0.5 case. The J p = 0.5 curves is shifted 
by T = —0.6 along the temperature axis to overlap with J p = 0.1 for comparison. The middle-value-based deviations 
S m are calculated for the gray area. The temperature resolution is AT = 0.0001. The width of gray window is 0.02. 


is derived in the same manner of standard derivation: 


8 m = 


\ 


1 N 

T T 'y X Xi ~ x i-m) 2 

V i=1 


( 21 ) 


The reliability of this deviation highly depends on the temperature resolution. Except the sharp transition range the 
results are acceptable with ST < 0.01 according to our experience. Fig[9] demonstrates the distributions of solutions 
with the J p = 0.1 and 0.5 in around T c range with resolution of AT = 0.0001. The J p = 0.5 curves is shifted by 
T = —0.6 for overlapping. The “middle-value-based” deviations of both cases are calculated in a selected range below 
T c . The results 5 m {J p = 0.5) = 0.001289 and S m (J p = 0.1) = 0.000419 agrees to the intuitive observation on the 
graph that higher J p corresponds to a larger fluctuation. 


D. Effect of 9 0 

So far the ground state angle 0 o is set to be 90°as a natural choice for the quadrilateral lattice. But it is possible to 
have situations with non-90° 0 Ol i.e. the quadrilateral cannot be the lowest energy state of particles. For example the 
six-carbon ring is known to be the most stable ring structure for organic carbon compound (e.g. benzene or cyclohex¬ 
ane), however the four-carbon ring is also thermodynamically possible as a metastable compound (e.g. cyclobutane 
or cyclobutadiene). We will investigate the effect of 9 a 90° on system’s thermodynamics in this section. 

The J p is fixed to be 0.2 for the 0 o variation. The 9 0 value of 30°, 60°, 90°, 120° and 150° have been investigated. 
The transition temperature and energy density at zero are summarized in table [II] The entropy behaviors are shown 
in FigfTOl 


TABLE II: The effects of 9 0 on T c and ground energy. 


(deg) T c E{T = 0) 

30 

3.622 -2.492 

60 

3.226 -2.255 

90 

3.096 -2.180 

120 

3.038 -2.147 

150 

3.011 -2.133 
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FIG. 10: The entropies with various 9 a = 30°, 60°, 90°, 120° and 150° with J p = 0.2. 
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e =30° 

o 

e o = 60 ° 
e =90° 

o 

0 = 120 ° 

o 

e =150° 



5 m (e o = 30°) = 0.002400 
S m (0 o = 60°) = 0.001771 
8 m (0 o = 90°) = 0.001060 
S m (0 o =12O°) = 0.000959 
8 m (0 o =15O°) = 0.000970 


FIG. 11: The solution distributions in around T c range of 9 0 = 30°, 60°, 90°, 120° and 150°. The shifts of the first 
four curves are T = —0.6, —0.21, —0.06 and —0.025 to overlap with the 150° for comparison. The middle-value-based 
deviations S m are calculated for the gray window of 0.005 width. The temperature resolution is AT = 0.0001. 


The results of 9 a = 90°, 120° and 150° basically agrees to the angle energy setup in Eq[2] and [3] that off-angle 
configurations will make the system less stable and lower the transition temperature. For 120° and 150° cases any 
randomized configuration must be an off-angle one therefore the lower T c and higher ground energy are expectable, 
although this effect is very small that in Fig[10]the three curves are very close. 

However 9 0 = 30° and 60° exhibit a higher T c and lower ground energy. That counter-intuitively implies if we 
force four particles, which is stable in equilateral triangle configuration, to form a quadrilateral, that will be more 
difficult to decompose. So far to the author’s knowledge there is no linear covalent bonds prefers an angle as small 
as 60°. Phosphorus white crystal has 60°bond angles but that is in 3D with regular tetrahedron structure. Due to 
lack of analogs in nature, it is not known if this phenomenon is realistic or it is simply a defect of our angle energy 
formulation. 

Fluctuation degrees of various 9 a are presented in the Fig[TT] calculated by Eq|2TI we have: 

S m (9 0 = 30°) = 0.002400, 

S m (9 0 = 60°) = 0.001771, 

S m (9 0 = 90°) = 0.001060, 

6 m (9 0 = 120°) = 0.000959, 

S m (9 0 = 150°) = 0.000970. 

The results follows the same principle in Fig[5]that the larger off-angle energy is, the larger the fluctuation is. 
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V. SUMMARY AND CONCLUSION 

We introduced the random angle variation in Husimi square lattice and generalized it to be a rhombus lattice, this 
is made possible by the independence property of recursive units. Considering a real system with coordination number 
of 4 in 2D, if we randomly sample four neighboring particles they are certainly in a quadrilateral configuration with 
random bond length and angles, even the system follows a deterministic dynamics. By this meaning we constructed 
a stochastic recursive lattice with Ising model on it, while no artificial randomness (the general “Chaos” or “Noise” 
in many models) on the dynamics is necessary to present the stochasticity in the system. 

For simplified treatment and easier calculation to explore the random lattice, the model is confined into one angle 
randomness and diagonal-proportional interaction. The Spontaneous Magnetization of ferromagnetic Ising model is 
presented in our model, but with a lower transition temperature T c since the off-angle configurations contribute a 
higher energy making the system less stable. The thermal fluctuation is induced by stochastic configurations and we 
introduced the quantity “middle-value-based” deviation S m to characterized the degree of fluctuation. 

The effects of the variation of diagonal interaction parameter J p and ground angle 9 a are investigated. Both T c and 
S m are behaved as expected based on the features of each variable, except the phenomena observed with 9 a < 90° are 
intriguing and may imply further interesting investigations. 
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